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Abstract. The algebraic classification of the Weyl tensor in arbitrary dimension n is 
recovered by means of the principal directions of its "superenergy" tensor. This point 
of view can be helpful in order to compute the Weyl aligned null directions explicitly, 
and permits to obtain the algebraic type of the Weyl tensor by computing the principal 
eigenvalue of rank-2 symmetric future tensors. The algebraic types compatible with 
states of intrinsic gravitational radiation can then be explored. The underlying ideas 
are general, so that a classification of arbitrary tensors in general dimension can be 
achieved. 



PACS numbers: 04.20.Cv, 04.50.-h, 02.40.Ky 
The Petrov classification (e.g. [321 1251 [5]) of 4- dimensional spacetimes can be 
reformulated by using the principal directions of the Bel- Robinson tensor. These are the 
causal vectors whose contraction with the Bel- Robinson tensor vanishes. The underlying 
ideas go back to [H El [121 E2] , are implicit in [25] and were fully exploited in [6] by using 
spinors. The result follows because the principal directions of the Bel-Robinson tensor 
coincide with the principal null directions of the Weyl tensor. To summarize, let [21191129] 

T a /3\^ = CapXaC^^ + C ap ^ a Ci3 p x a — -g a /3gx^Cp T(TU C pTal ' , if n = 4 (1) 

be the Bel- Robinson tensor of a 4- dimensional spacetime, where C ap \ a is the Weyl tensor 
and gxfj, the metric tensor. T a /3x^ is completely symmetric and traceless. The Petrov 
classification can be described as follows [6]: there exists a null vector l p such that 

(i) T apx ^ = Petrov type N 

(ii) T a px^ X ^ = but T afi x^ ± Petrov type III 

(hi) T al 3x^e X ^ = but T al 3x^ X ^ 1 ^ Petrov type II or D 

(iv) Tapx^iPi^ = but T a /3x^£ x ^ Petrov type I 

Furthermore, one can distinguish types II and D by using 

(v) there exist two linearly independent null vectors £ M and such that T a px^^ X ^ = 
and T a px^k x k p = Petrov type D . 
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Actually, one can drop the assumption that the vectors are null by just assuming that 
they are causal (i.e., timelike or null). 

In recent years, the algebraic classification of the Weyl tensor in higher dimensions 
has been constructed based on the existence of WANDs (Weyl aligned null directions) 
and their alignment order (TTJ [23], see pU] for a review. Criteria to characterize 
the (multiple) WANDs, analogous to the classical Bel-Debever criteria [U |5l [131 Eh 
have been also recently obtained by Ortaggio [21]. In 4-dimensonal spacetimes, these 
criteria are known to be equivalent to the number of times one has to contract with 
the Bel-Robinson tensor, and thereby they provide an alternative view of the Petrov 
classification. I am going to show that the same happens in arbitrary dimension. 

To that end, one only needs the correct generalization of the Bel-Robinson tensor 
to higher dimensions. As argued in [29], this is the so-called 'superenergy' tensor T{C} 
of the Weyl tensoi[||, whose explicit expression in arbitrary dimension n is 



1 



TafiXfi = C a p\o-Cf3 P n + C ap^o-C [3 P \ — —g a {jC p T \ a C P 



/' 



1 1 

- -g\fiC apcTT Cf3 pcrT + -g a /3gx^Cp Tau C pTcrl/ (2) 

2 o 
from where one deduces the symmetry properties 

T a p\p = T( a (3)(\n) = Txpa/3. (3) 

However, T a p\p is completely symmetric only in dimensions n 6 {4,5} [29] . Similarly, 
unlike in the case n — 4, all traces of the tensor (J2J) are (generically) different from zero 
for all n > 4. Expression ([2]) reduces to the original expression ([1]) if n = 4 due to a 
4-dimensional identity [29| [T6]. 

As any other superenergy tensor, the generalized Bel- Robinson tensor (T2]) is a future 
tensor, that is to say, it satisfies the dominant property [29] 

T^aX^sM > 0, (4) 

for arbitrary future-pointing vectors v", v%, v$ and v p . Inequality @ is strict if , v% , v£ 
and v% are all timelike [291 [7]. Thus, any possible causal vector v p with the property 

7^ v v Vv V = (5) 

must be null. Causal vectors satisfying (jSJ) define the principal directions of Ta/3\p, 
[3 UHl [20]. As shown in [26] with full generality -for the superenergy tensor T{A} 
of any tensor A and in arbitrary dimension — , these are precisely the principal null 
directions of the Weyl tensor, that is to say, those null vectors satisfying 

V[pC a]pa[x Vp]V p v a = 0. 

This is, in fact, the characterization of WANDs [23| 127] . 

J Given any tensor A the super-energy construction [55] is a general method to build a — essentially 
unique — tensor TjA} quadratic in A and future. Future tensors are those satisfying the dominant 
property (e.g. Q); see [29] section 4, section 2 in [7], or the Appendix in [20] for further details. Due 
to historical reasons [lj [2 [H SI El [12l [13l H31 [22l [29], T{A} is called the super-energy tensor of A. 
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Proposition 1. At any point of a causally orientablen- dimensional Lorentzian manifold 
where the Weyl tensor does not vanish, a causal vector £ p satisfies: 

(i) Tapx^tHHi* = if and only if £[pC a]pa[ x£ p] £ p £ a = , that is V is a WAND, 
(n) T a px^ a in x = if and only if t^C a]p[aX £ ^ = 0. 
(Hi) %p X ^ = if and only if £ [p C a]pa x£ p = 0. 

(iv) T a /3\tj,£ a £ x = if and only if £[/3C ap \[ a x£ p ] = and C ap [ a \t^ p = (which together 

actually entail £[pC a ] pa x£ p — 0j. 
(v) T a /3\^ a = if and only if £^C ap]aX = 0. 
In all cases £^ is necessarily null. 

Proof. As stated above, the first point (i) is a consequence of a fully general result for 
superenergy tensors T{A} [26]: once one has sorted out the number of antisymmetric 
blocks of the seed tensor A, the principal directions of are precisely the null 

directions such that their contraction (or inner product) followed by the exterior product 
on all skew-symmetric blocks of A vanishes. However, I am going to include an 
elementary proof here for the case of the Weyl tensor as this will be illustrative and 
because it will be useful in the rest of the cases. 

Take any null £ p '. Contracting with (j2J) one gets 

T^F&lH* = 2C pa C p ° 

where I have defined 

Cap{£) = c ap pjn°. (6) 

Taking into account the obvious properties C a p = Cp a and C a p$ = it follows that 
C a pC a ^ = is equivalent to (for instance by decomposing on a basis) 

C a(5 = l a v p + e p v a (7) 

for some vp (such that £ p v p = 0). But this means that £[xCd\\p^^] — 0, proving (i). 

To prove (ii), suppose that T a pxp,£ a £^£ X = 0, so that in particular (J7|) must hold, 
due to (i). Then a direct computation using (JHJ) and (I7j) repeatedly gives 

= T aP x,£ a £ p £ x = -i* \2v p v p + ho pra C^j (8) 
where I have set 

Cpxn(£) = C a /3x^£ a 

with the obvious properties Cp Xp = C^j, C [l3Xp ] = 0, £ p Cp Xp = 0, £ x Cp Xli = 
Cpp = £pv p + £fjVp. Choose another null vector such that £ fl k^ = 1, define 
F Xp = k p C p x p = F[\n] an d put a hat on any tensor orthogonal to the timelike plane 
spanned by £^ and k p . Then a typical decomposition proves 

CW = + hAp p - £ p A/3x + Cpx» (9) 
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with lf>A pu = kPA pa = kPA ap = 0, £PA ap = -v a + (kPv p )£ a and £PF pu = v a + (kPv p )£ a . 
A direct calculation shows that ([8]) becomes 

-v P v p - ^C^C^j = 

which implies C pra = and = \A£p \ Then, (0) simplifies to C pu = A£ p £ u and ([9]) to 
= £fsF Xll + £ x Ap p - £ p Ap X (10) 
from where £\ a Cp\\xp,£v\ = 0, proving (ii). Observe that C[pxn] = implies now 

F^v = 2A [pu] where = £ p p v - £ v p p + F pv (11) 

for some pp. 

Assume now T a f3x^£ a £^ = so that in particular (fit)]) and ( TTTj) hold. Using (|T0|) a 
direct computation gives 

o = T aP x/ a £ p = £x£,A pa Ap° 

so that A pa = hence C pXf i = £pF Xp (with F Xp = 0). Thus £[ a C^ Xpi = proving (iii). 

If on the other hand T a /3x P £ a £ x = 0, and using repeatedly (ITUj) and (ITT]) which 
remain valid in this situation, the calculation provides 

= T apx »£ a £ x = tplp (f p Jp° - A pa A°P + ic^C^J (12) 

so one needs to compute the last term. At this stage, by setting qp = (1/2) (£ p p p )kP — pp 
the Weyl tensor takes the form 

C af 3 Xf i = (£ a kp - £pk a )(£ x q p - £ p q x ) + (£ a qp ~ £pq a )(£\k p - £ p k x ) 
-(£ a kp - £f}k a )Fx p -F a p(£xk p - £ p k x ) + Ak[ X A p ] [l3 £ a] + Ak[ a A^ p £ x \ 

+2£[aUj3]X l i + 2£[\U p ]al3 + ^[aYpW^X] + C a px P 

for some Upx^ — ^s[A/i], Vp p = V p p and Capx^ with the same symmetry properties as the 
Weyl tensor. It follows that 

C a( 3x,C a ? x P = A(£ p pPf - AF pa FP° + sA^p+C^C^p 

so that (j!2p leads to 

A(£ pP Pf + AF^FP* + C^C*^ = 

implying £ p p p = 0, = and Capx/j, — 0. Using the trace-free property of the Weyl 
tensor, this last result implies also that A( Xp ) = 0, which together with Fx p = provides 
A Xp = 0. All this means £[pC ap \\ a x£p\ = (and furthermore C pu = 0, as A = —£ p p p = 0) 
plus the condition found previously in (iii), proving (iv). Observe that (iv) always 
implies (iii). 



Algebraic classification based on the superenergy tensor 



5 



Finally, assume that T a p\^i a = so that everything that has been derived in (hi) 
and (iv) also holds, that is, A^ v = 0, F^ u = 0, £ p p p = (so that p p = p p ) and C a/ 3x p = 0. 
One easily gets 

= T a ^£ a = -£p£ x £» (p p p p + \u P r,u pTU ^j 

providing p p = and Upxn = so that the Weyl tensor adopts the very simple form 
C a p\fi = ^[aV/3][^x\ and thus £[ v C a ^x P = which ends the proof. ■ 

Combining this proposition with the results in [21] one immediately obtains the 
algebraic classification of the Weyl tensor using the notation of [HI [10]. If C a p XfJ ,\ x ^ 0, 
then its algebraic type can be characterized according to whether there exists a null 
vector £ p at x such that: 

N -<=>- Tapx^ = 0. In this case, £ p is the unique null vector (up to proportionality 
factors) with this property, and it defines the unique WAND. 

Ill T a px^£ a £ x = but T a p Xfl £ a ^ 0. Observe that in this case T a px^£ a £ p = holds 
necessarily. The null £ p with this property is uniquely defined and characterizes the 
unique multiple WAND. 

II or D T a p X / a ^? X = but T aP x/ a £ x is non-zero for all such £ p . 

I <^ Tapx^e 13 ^ = but T a px^ a £ P £ X + for all such £". 

G <^=> T a px„£ a £ P £ X £^ ^ for all null £". 

In addition, one can distinguish between types D and II using the following 
D there exist two linearly independent null vectors £ p and k p satisfying 

T a pxp,£ a £^ £ X = %px^k a k^k x = (observe that then T a px/ a £ X and T a/ 3x^k a k x must 
be different from zero). 

For these cases D and II, one can further characterize some of their subcases as follows. 
If T a px^£ a £^ = 0, then the types are ll a bd or D a bd, in the latter case T a pxpk a k^ = 
actually holds too. Note that the cases with T a /3x^£ c = and T a i3x^k a k x = are 
impossible. 

Finally, for types I, II (Il a bd) and III there exist secondary subtypes Ij, Ilj (Iliabd) 
and Illj which are simply characterized by the existence of another simple WAND k p , 
that is, a null vector k p which is linearly independent of the given £ p and such that 
T a ^k a k^k x k p = 0. 

One can now explore the algebraic types compatible with a state of intrinsic 
gravitational radiation at any point x. Following Bel [U [5] this will happen whenever 
T a pXfj, uau u is a non-zero null vector for all timelike vector u p . This definition seems 
satisfactory in arbitrary dimension because in static cases, when there exists a timelike 
hyper surf ace-orthogonal Killing vector there will never be intrinsic gravitational 
radiation. To see this, use the result in [21] stating that for such T a px^ a ^£ x = 
(in Ricci-flat cases for simplicity). This rules out types G, Ij and D according to [28] . 
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Apart from its intrinsic interest, the previous characterization of the algebraic types 
of the Weyl tensor can be useful in order to compute the WANDs and to actually classify 
explicit spacetimes. To take full advantage of this alternative one must use some of the 
general properties of future tensors, see [H [20] for the needed details. 

Using that T a p xn is a future tensor, first of all one can compute T p px l J' p aT u- If the 
result is zero, then necessarily T a px P = t a tpxp, for a null £^ |7| and thus £ a T a p\p = 0, 
so that the Petrov type is N. Observe that one does not need to know anything about 
WANDs to check this result. In order to know the unique WAND, simply contract with 
an arbitrary timelike vector w M thrice: = T a ^u a u^v? ' . 

Suppose then that TppxpCf 9 cxtv does not vanish. For any causal (timelike or null) w^ 1 
define Tp^w) = T a pxj_iUJ a w x which is a symmetric rank-2 future tensor. If Tpp{w) = 
then u> M is necessarily null and the Petrov type is III. If Tp^w) ^ 0, then one looks 
for its null eigenvectors. In order to know them only the principal eigenvalue, call it 
A, is needed [7J [20] . One can calculate A by classical methods — by power iteration 
and/or solving an algebraic equation. If A is degenerate or corresponds to a double null 
eigenvector then [7J [20] Tp^i 13 = A£ /t for some null i^, so that Tp^i^ = and one 
derives T a px^^ = leading again to type III. The same procedure can be performed 
with the symmetric rank-2 future tensor Tx^iu) = T a pxpU a u^ for an arbitrary timelike 
u*. If T x / X = X£p for some null V then Tx/ X l» = and one obtains T a px^ x ^ = 0. All 
null vectors with this property are those in the eigenspace associated to the principal 
eigenvalue A of Tx^u). If there are at least two of them linearly independent, then the 
Petrov type is D a bd- If there is only one, the Petrov type is ll a bd- Observe that, for the 
iteration procedure to actually finding these princiapl null eigenspaces, one can use as 
many Txn(u), with different u^ 1 , as desired. Actually, one can even use tensors of type 
T a i3Xfj,u a v 13 for any two timelike and v 1 *. 

If neither Txfj,(u) nor Tp^u) have null eigendirections — which happens when A 
corresponds to a non-degenerate eigenvalue with timelike eigenvector — then one knows 
that T a px^ X ^ ^ and T a px^^ ^ for all possible null I*. Choose then an 
arbitrary null £ M and construct Tx^i) = T a px^ a ^, which is a non- vanishing future 
tensor. Computing its principal eigenvalue — which depends on — , if it happens to 
be non-degenerate with timelike eigendirection, then there are no WANDs. If it is either 
degenerate or corresponds to a double null eigenvector but does not (respectively does) 
vanish, then one must check if the corresponding null eigenvector coincides with l 9 for 
some choice of the latter, in which case T a pxfi IH^lH* = (resp. T a px^ a ^^ = 0). This 
takes care of types D, II, I and G. The calculations in these situations may be long. 

It is also possible to provide characterizations of the different types by using tensors 
and scalars obtained by taking powers of T a pxp, and then contracting some or all of the 
indices. In 4-dimensional spacetimes these results are known but not easy to derive, see 
[El El [T7], HE] • In higher dimensions they may be even more involved. Nevertheless, they 
would be very important providing invariant ways of determining the Weyl algebraic 
types. In a similar vein, the invariants and concomitants written in terms of a general 
electric-magnetic decomposition could be used, see [H El [5j El El EQ] . 
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Notice that the method outlined here applies, mutatis mutandis, to the Riemann 
tensor, and actually to any double (2,2)-form without the symmetry between pairs of 
indices. Actually, I would like to remark that an algebraic classification of any tensor 
A can be achieved by the same method, on using its superenergy tensor T{A}. In 
general this is a tensor with 2r indices, distributed in r symmetric pairs according to 
the number r of anty-symmetric blocks of indices of the seed tensor A [291 CSj- The 
number of different algebraic types depends on r. As a sufficient illustrative example, 
consider the case of a rank-2 tensor A^ u (no symmetries assumed), whose superenergy 
tensor is given by [29] 



with the property T a(3X ^ = T ia(3){Xfl) ( if A fw happens to be symmetric then T a p Xfl = T Xmfj 
holds too. In the skew-symmetric case the proper superenergy tensor has only a pair 
of indices, and the tensor above is a generalization containing it [31]). An elementary 
calculation provides the following equivalences (in all cases is null): 



• T a p Xll i a inH» = o A^ u m v = 

• T^PlH* = e»A^j T] = 

• T^lHH* = £ [tl A T]v t = 

• TW*^ = ^ PA^ = 

• 7W A *" = ^ A^t = 



As a final remark, I would like to mention that these classifications can be 
refined, obtaining more information involving several null directions simultaneously, 
by considering a generalization called the 'mathematical energy tensor' [3T], but this is 
out of the scope of this short communication. 
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